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Abstract — In this paper, the model of the Photovoltaic (PV) 
array has been explained in which the basic current and voltage 
equations for a PV cell have been defined. It is discussed that 
output response of PV system depends upon incident solar 
radiation and ambient temperature. The modeling of solar PV 
single stage grid connected system at different power factors is 
done. The nature of real power generated by solar PV array 
through converter has been shown and proved that, whenever 
the power from grid is un available, the real power requirement 
of the load is achieved by converter. Maximum Power point 
Tracking (MPPT) technique in which data through look up 
table is used through which behavior of actual DC link voltage is 
discussed. Finally, the simulation of proposed system is 
performed under two different cases. In first case, the effect of 
changing power factor on active power, reactive power and 
Total Harmonic Distortion (THD) values is observed. It is found 
that the THD of grid current increases with increase in the 
phase angle of grid current converter voltage. It affects the 
active power flow among converter, load and utility grid. In the 
second case, the effect of changing frequency on active power, 
reactive power and THD values is noticed. It is observed that the 
THD of converter current increases, whereas, the THD of grid 
current remains constant. 

Index Terms — Photovoltaic, harmonic, power factor, power 
and frequency 


I. INTRODUCTION 

The non-renewable energy sources and fossil fuels 
such as coal, petroleum and natural gas have been extensively 
employed for thousands of years on the earth. After their 
burning, these fossil fuel produces excessive heat which 
results severe environmental pollutions such as air pollution, 
water pollution, and soil pollution. The fast depletion of the 
fossil fuels may lead to shortage of ever increasing electrical 
power demand of the world for various domestic and 
industrial applications. This leads electric power engineers 
and researchers to think for alternative available renewable 
resources such as wind, Photovoltaic (PV), Fuel Cells (FC), 
etc. 

The conventional Maximum Power Point Tracking 
(MPPT) methods can be classified into number of categories 
as mentioned in [1]. From the literature survey, the first 
category consists of constant voltage or current reference 
methods. In this, the controller sets the voltage or current 
reference to an empirical ratio of the open circuit voltage or 
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short circuit current. Although it is very simple and quick 
method, but can’t extract maximum power under all 
environmental conditions and has a low tracking efficiency 

[2] . The second category consists of the Perturbation and 
Observation (P&O) method. P&O MPPT method use the idea 
of hill climbing algorithm. It works well, when the solar 
radiation does not vary rapidly with time. Due to simplicity 
and ease of implementation, it is most widely applied methods 
in PV industry. However, their drawbacks are slow tracking 
speed and oscillations around Maximum Power Point (MPP) 

[3] . To reduce or eliminate these drawbacks, P&O methods 
have been improved by using variable step-sizes [4]. 

Another disadvantage with P&O MPPT method is 
that when the environmental conditions change rapidly, i.e. 
when the solar radiation changes rapidly, the operating point 
may deviate from the MPP by tracking towards a wrong 
direction. This is because the algorithms are not able to 
distinguish between changes caused by the perturbation or by 
the changing environmental conditions. Various methods to 
overcome this problem have been proposed in references 
[5-6]. The third MPPT category is the Incremental 
Conductance (IC) method [3]. This method tracks faster and 
has a better performance than the P&O method, especially 
under rapidly changing atmospheric conditions. There are 
also some other MPPT algorithms. For example, an MPPT 
method which senses the output voltage or current of 
converter is presented in [7]. The Pulse Width Modulation 
(PWM) based converter is modeled as a time-varied 
transformer, so that the input voltage or current can be 
calculated from the output voltage or current. By analyzing 
linear and non-linear loads, it was found that for loads with 
non-negative incremental impedance, the output current and 
voltage increases monotonically as the output power 
increases. Thus maximum output current (or voltage) implies 
maximum output power. Therefore, the MPPT algorithm can 
be implemented based only on a single output parameter, 
which can simplify hardware design and control algorithm. 
With this method, it is the output power of Voltage Source 
Converter (VSC) that is maximized, rather than the solar PV 
array power. Thus solar PV array may be operating slightly 
off the MPP, but the total efficiency is optimized. The 
characteristics of solar radiation are found to be constantly 
variable. The atmospheric conditions, climate, and 
geographic characteristics are most important parameters, 
which determine the solar radiation quantity received on a 
given point of the earth [8-10]. 

II. BASIC MODELING OF SOLAR PV SYSTEM 

As mentioned in reference [11], for a known 
temperature and a known solar radiation level, a model is 
obtained which is then modified to handle different cases of 
temperature and solar radiation levels. The solar cell 
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operating temperature varies as a function of solar radiation 
level, and ambient temperature. The variable ambient 
temperature affects the cell output voltage, and cell 
photocurrent. These effects are represented in the model by 
Equations (l)-(2) as: 

C n =l + /3,(T a -T x ) (1) 

(2) 

A 

where, p t = 0.0042 and y = 0.062 for the cell used, and T a = 



(a) 


200 °C, is ambient temperature during the cell testing. This is 
used to obtain the modified model of cell for another ambient 
temperature T x . A change in solar radiation causes a change in 
cell photocurrent, and operating temperature, which in turn 
affects cell output voltage. If solar radiation increases from 
S xl to S x2 , cell operating temperature and photocurrent will 
also increase from T xl to T x2 and from I ph] to I ph2 , respectively. 
Thus change in operating temperature and photocurrent due 
to variation in solar radiation are expressed as, 

C„=l + /3,a,(S x -S c ) O) 


where, S c is benchmark reference solar radiation level during 
cell testing to obtain modified cell model. S x is new level of 
solar radiation. The constant a s represents slope of change in 
cell operating temperature due to change in solar radiation 
level and is equal to 0.011 for solar cells used. Using 
Equations (l)-(4), the new values of cell output voltage V an 
and photocurrent I phn at new temperature T x and solar 
radiation S x are obtained as, 

Van = Ctv C sv V a (5) 

Iphn ~ Cfi C si Ip r; (6) 

V a and I ph are benchmark reference cell output voltage and 
reference cell photocurrent, respectively. 

The I-V and P-V curves for changing solar radiation 
but a fixed temperature (25 °C) is shown in Figure 1 (a) and 
Figure 1(b), respectively. The I-V curve it is found that there 
are two regions in the curve: one is current source region 
whereas other is voltage source region. In the voltage source 
region (in right side of the curve), the internal impedance is 
low whereas, in current source region (in the left side of the 
curve), the impedance is high. Solar radiation and ambient 
temperature plays an important role in predicting the I-V 
characteristic. The effects of both factors are considered 
during the designing of PV system. Whereas the solar 
radiation affects the output, temperature mainly affects the 
terminal voltage. Figure 2 (a) and Figure 2 (b) gives the 
simulated I-V and P-V characteristic for various temperatures 
at a fixed irradiance at 1000 W/m 2 , respectively. 



VOLTAGE <V) 


(b) 

Figure 1 (a) I-V and (b) P-V characteristic, of a solar 
array for a fixed temperature but varying irradiance 




(b) 

Figure 2 (a) I-V and (b) P-V characteristic of a PV array 
under a fixed irradiance but varying temperatures 

A. CONFIGURATION OF SOLAR PV SYSTEM 
As shown in Figure 3, the Power Conditioning System 
(PCS) is comprised of VSC along with its control, MPPT 
control; Phase Focked Foop (PFF) and LCL filter for the 
removal of harmonics. PFF algorithm is used for grid 
synchronization. The output of the VSC is connected to grid 
through a inductance so that the flow of power can take place 
on both sides. The presence of inductance creates a phase 
angle which is responsible for the flow of active and reactive 
from both sides. A look-up table is used in Matlab simulation 
in which the data is fed for the voltage at MPP. The data 
which has been used is 600, 700, 800, 900, 1000,1100 and 
1200 V at various solar radiation range 88, 90, 92, 94, 96, 98 
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and 100 W/m 2 . No. of series and parallel connected PV cells 
are 2200 and 26, respectively. The three phase series RLC 
load used is 440 V,50 Hz, P=63 kW, Q=20 kW (+ve). DC 
voltage regulator gains are 0.45 and 18, whereas DC current 
regulator gains are 0.026 and 5. 



Figure 3 Simulated solar PV system connected with load 


III. SIMULATION RESULTS AND DISCUSSION 

For the proposed system, the total simulation time 
period is 0.3 s. Figure 4 (a) shows the nature of the solar 
radiation which acts as input to the proposed system. Since 
array current is directly proportional to the solar radiation, its 
nature is approximately proportional to the solar radiation 
curve, Figure 4 (b). Figure 4 (c) gives the PV array voltage, 
whereas Figure 4 (d) shows the PV array power. The output 
power of PV array is product of PV array voltage and current. 
It is possible to observe that the approximate value of the PV 
array power is 48 kW. 




(b) 




(d) 


Figure 4 (a) Solar radiation (b) Solar array current (c) 
Solar array voltage (d) Solar array power 

Figure 5 (a) depicts the nature of actual DC link voltage 
around the reference MPPT voltage. The blue line indicates 
the tracking of actual DC link voltage around reference 
voltage. As the ambient temperature T x increases and solar 
radiation increases, the actual value closely follows the 
reference voltage. The initial transients can be controlled by 
the further tuning of PI controllers. Figure 5 (b) shows that as 
the solar radiation level increases the grid power decreases. 
The power requirement of the load is met by PV array 
through VSC converter. This has been shown by Figure 5 
(c), which depicts the VSC power. Figure 5 (d) shows the 
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(d) 

Figure 5 (a) Actual DC link voltage and MPPT 
reference voltage (b) Grid output power (c) Inverter 
output power (d Modulation index 

The performance of the system is evaluated at 
different power factors and frequencies. The total simulation 
period has been chosen as t = 0.35 s. The following are the 
two case studies: 

Case I at different power factors and same frequency; Case II 
at different frequencies and same power factors. 

A. Case I at different power factors and same frequency 
•at unity power factor , 440 V, 50 Hz 

As shown in Figure 6 (a), the active power is 
generated at different temperature and solar radiation, at 
unity power factor. The system has been tested at constant 
frequency of 50 Hz. This active power has been generated 
according to data generated through MPPT used through 
look-up table. The grid power completes the load power 
demand. As it decreases, the load power demand is met with 
the solar PV array through the IGBT based converter system. 
The initial transient at t = 0.02 s can be more controlled by 
introducing tuning of PI controllers. Similarly, the reactive 
power response has been shown in Figure 6 (b). Reactive 
power compensation of the load is done by solar PV through 
converter system (shown by blue line). At t = 0.05 s, the 
reactive power generated by grid decreases till t = 0.35 s. 
The reactive power of the load is shown as constant at 20 
kVAR. 
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THD values from Fast Fourier Transform (FFT) analysis for 
VSC current and grid injected current are found to be 1.05 % 
and 0.49 %, respectively. 

• at power factor = 0.866, 440 V, 50 Hz 

As depicted by Figure 7 (a), there has been a non-uniformity 
in the active power generation between the period t = 0.07 s to 
t - 0.1 s, which is due to decrease in the power factor from 
unity to 0.866. This is due to the fact, that the power factor is 
directly proportional to active power. Similar behavior has 
been shown by reactive power response between the period t 
= 0.02 s to t - 0.1 s, Figure 7 (b). 
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(b) 

Figure 7 (a) Active power (b) Reactive power, variation of 
VSC, load and utility grid (Power factor=0.866, 

50 Hz) 

The THD response of VSC current and grid injected 
current at 0.866 power factor, respectively. The THD values 
are found to be 7.19 % and 2.71 %, respectively. Thus there 
has been the increase in THD with decrease in power factor. 

•at power factor = 0.707, 440 V, 50 Hz 

At 0.707 power factor, the response becomes more 
oscillatory as shown by active power generated by solar PV 
through converter system, Figure 8 (a). The maximum power 
generated by utility grid is 80 KW. However, there has been 
decrease in the overshoots in the reactive power response 
after t = 0.15 s, Figure 8 (b). 


(b) 

Figure 6 (a) Active power (b) Reactive power, variation of 
VSC, load and utility grid (Unity power factor, 50 Hz) 


utility grid power 


(a) 
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(b) 

Figure 8 (a) Active power (b) Reactive power, variation of 
VSC, load and utility grid (Power factor=0.707, 

50 Hz) 

The THD values are found to be 8.28 % and 2.67 %, 
respectively. Again, there has been the increase in the THD of 
VSC current with decrease in power factor. 

• at power factor = 0.5, 440 V, 50 Hz 
Figure 9 (a) depicts that the overshoot increases to 1.2 x 10 5 . 
Thus the oscillatory response increases with the decrease in 
operating power factor. From Figure 9 (b), the reactive power 
generated by utility grid is 1.1 x 10 5 . However, in this case, 
the compensation is met with the grid power from 0.04 s to 
0.08 s, thereafter, completed by reactive power generated by 
VSC. 



(a) 



(b) 

Figure 9 (a) Active power (b) Reactive power, variation of 
VSC, load and utility grid (Power factor=0.50, 50 Hz) 


The THD values are found to be 9.68% and 2.58%, 
respectively. Thus, it has been found that there has been the 
increase in the THD of VSC current from 1.05% to 9.68% as 
the power factor decreases. 


B. Case II at different frequency and same power factors 

This sub-section presents the behavior of proposed system 
at unity power factor and two different frequencies. 

• at unity power factor, 440 V, 49 Hz 

Figure 10 (a) depicts the active power response 
whereas; Figure 10 (b) depicts the reactive power response of 
VSC, load and utility grid. The THD values of VSC current 
and grid current are found to be 0.49% and 1.84%, 
respectively. 



(b) 

Figure 10 (a) Active power (b) Reactive power, variation 
of VSC, load and utility grid (Unity power factor, 49 Hz) 


• at unity power factor, 440 V, 51 Hz 

Figure 11 (a) depicts the active power response 
whereas; Figure 11 (b) depicts the reactive power response 
of VSC, load and utility grid. The THD values of VSC 
current and grid current are found to be 6.52% and 1.84%, 
respectively. This it has is concluded that change in 
frequency has no effect on the THD of grid current whereas, 
THD of VSC current has been increased from 0.49% to 
6.52%. 



(a) 
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(b) 

Figure 11 (a) Active power (b) Reactive power, variation 
of VSC, load and utility grid (Unity power factor, 51 Hz) 


IV. CONCLUSION 

In this paper, the simulation of the proposed system is 
performed under two different cases. In first case, the effect of 
changing power factor on active power, reactive power and 
THD values is observed. It is observed that the THD of grid 
current increases with increase in the phase angle of grid 
current VSC voltage. It affects the active power flow among 
VSC, load and utility grid. In the second case, the effect of 
changing frequency on active power, reactive power and 
THD values is noticed. It is observed that the THD of VSC 
current increases whereas, the THD of grid current remains 
constant. 
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